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ABSTRACT 


"  the  mechanical  properties  of  3  types  of  cellular  solids  (flexible, 
plastic  and  brittle)  have  been  measured  as  a  function  of  density.  The 
results  are  compared  with  models  for  the  stiffness,  strength  and  den- 
sification;  and  constitutive  laws  are  developed.  Data  and  models  for 
each  type  of  cellular  solid  are  combined  to  develop  mechanism-mode 
maps  which  summarise  the  properties  in  a  single  diagram.  Natural 
cellular  materials  fit  the  same  pattern;  maps  are  presented,  as  an 
example,  for  wood.  The  maps  help  in  design  and  in  the  selection  of 
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TABLE  1:  SYMBOLS  AND  UNITS 


a 

e 

P 

P  (e) 
ps 

p/p5 

E* 

E 

s 

* 


t 

l 


i  Ce) 


Applied  compressive  stress  Of/m2) 

Nominal  compressive  strain  (-) 

Initial  density  of  cellular  solid  (kg/m3) 

Density  after  compressive  strain  e  (kg/m3) 

Density  of  cell  wall  material  (kg/m3) 

Relative  density  (-) 

Young’s  modulus  of  cellular  solid  (N/m2) 

Young's  modulus  of  cell-wall  material  (N/m2) 


Elastic  collapse  or  plateau  stress  of  elastomeric  (N/m2) 
foam 

Plastic  collapse  or  plateau  stress  of  plastic  foam  (N/m2) 


Yield  strength  of  cell-wall  material  (N/m2) 

Crushing  stress  of  brittle  foam  (N/m2) 

Modulus  of  rupture  of  cell-wall  material  (N/m2) 

Cell  wall  thickness  (m) 

Cell  size  or  cell  wall  length  (m) 

Length  of  uncollapsed  members  at  strain  c  (m) 


iQ  Initial  length  of  an  uncollapsed  member  (m) 

I  Second  moment  of  area  of  cell  wall  or  edge  (m* ) 

F  Force  acting  on  a  cell  wall  (N) 

F  Euler  buckling  load  for  cell  wall  (N) 

cr 

M  Fully  plastic  moment  of  cell  wall  (Nm) 

P 

M^  Moment  which  will  just  fracture  cell  wall  (Nm) 

Ci  -  Ce  Dimensionless  constants  (~) 

h  Initial  height  of  sample  (m) 

o 

h  Height  after  strain  e  (m) 

P  Load  00 


1. 


INTRODUCTION 


Polymeric  foams  have  certain  characteristic  mechanical  properties. 
Elastomeric  or  flexible  foams,  in  simple  compression,  are  linear-elastic 
to  a  strain  of  about  5  Z.  Then  the  cell  walls  buckle  and  the  foam  col¬ 
lapses  at  a  nearly  constant  stress  (giving  a  non-linear  elastic  deforma¬ 
tion)  until  the  cell  walls  touch  and  the  stress-strain  curve  rises  steeply 
Rigid  polymers  and  metals,  when  foamed,  have  a  similar  stress-strain  curve 
but  for  a  different  reason.  Like  flexible  foams,  they  are  linear-elastic 
to  a  strain  of  roughly  5  %.  Then  they  suffer  plastic  collapse,  compres¬ 
sing  plastically  at  a  nearly  constant  stress  until  the  cell  walls  touch, 
and  the  stress-strain  curve  rises  steeply.  Brittle  foams,  too,  show  an 
initial  linear-elastic  regime.  But  when  the  stress  is  reached  at  which 
cell  walls  fracture,  the  stress-strain  curve  becomes  irregular  (though 
roughly  horizontal)  and  the  foam  crushes  at  (roughly)  constant  load. 

Fig.  1  shows,  schematically,  the  structure  of  cellular  solids. 

Some  have  open  cells :  the  solid  material  is  distributed  as  little 
beams  which  form  the  cell  edges.  Others  have  closed  cells',  the  solid 
is  distributed  as  little  plates  which  form  the  cell  faces.  The  mechan¬ 
ical  properties  reflect,  to  some  extent,  this  distribution.  In  practice, 
most  man-made  foams  (even  those  with  closed  cell  faces)  behave  as  if  they 
had  open  cells  because  surface  tension  draws  much  of  the  solid  material 
into  the  cell-edges  during  manufacture.  For  this  reason,  we  discuss  open¬ 
cell  foams  in  detail,  but  treat  closed  cell  foams  only  in  passing. 

The  mechanical  properties  of  a  cellular  solid  can  be  related  to 
the  mechanics  of  bending,  buckling,  plastic  collapse  and  brittle  frac¬ 
ture  of  its  cell  walls.  Each  part  of  the  stress-strain  curve  can  be 
modelled  (1,2, 3, 4).  The  models  give  equations  for  Young's  modulus  E*» 


the  elastic  buckling  stress  o^,  the  plastic  collapse  stress  a^,  the 
crushing  strength  a^,  and  terminal  rise  in  strength,  in  terms  of  the 
density  and  properties  of  the  material  of  which  the  foam  is  made. 

This  paper  seeks  to  test  and  extend  these  models,  (checking  them  against 
data  from  the  literature  and  new  data,  described  below),  to  derive  constitutive 
laws  for  design  with  cellular  solids,  and  to  develop  diagrams  which  summarise 
the  overall  mechanical  response  of  each  type.  Symbols  are  defined  in  Table  1. 


2.  EXPERIMENTAL  RESULTS 


2.1  Foamed  Plastics  and  Ceramics 


We  tested  samples  of  commercial  flexible  foams  (a  polyethylene  and  a 
polyurethane) ,  samples  of  a  commercial  rigid  foam  (a  polymethylacrylimid) 
and  an  experimental  batch  of  a  brittle  ceramic  foam  (mullite).  The  materials 
their  origins  and  the  properties  of  the  cell  walls  are  listed  in  Table  2. 
Their  structures  are  shown  in  Fig.  2. 

Compression  tests  were  carried  out  on  blocks  of  foam  of  a  convenient 

o 

size  (the  size  depending  on  the  stiffness) ,  at  a  temperature  of  18  C  and 
a  strain-rate  of  roughly  2  *  10~Vs.  Results  are  shown  in  Figs.  3  to  6. 

The  axes  are  the  nominal  stress  (the  load  F  divided  by  the  initial 


section  A  ) : 

o 


and  the  nominal  compressive  strain: 
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where  h  is  the  height  of  the  sample  after  a  strain  e,  and  hQ  is 
the  original  height.  When  foams  are  compressed  beyond  a  strain  of  a 
few  %,  there  is  almost  no  lateral  spreading,  so  the  nominal  and  true 
stresses  are,  for  all  practical  purposes,  identical.  The  nominal 
compressive  strain  is,  of  course,  limited  to  the  range  0  to  1. 

All  the  stress-strain  curves  show  three  regions:  a  linear  elastic 
region ;  a  long  plateau  where  the  stress  is  almost  independent  of  strain; 
and  (for  all  but  the  brittle  foams)  a  final  region  of  densification  in 
which  the  stress-strain  curve  rises  steeply.  Young's  modulus  E  *  of 
the  rigid  foams  was  meaured  by  using  clip  gauges.  The  density  p  of 
each  foam  was  measured  by  conventional  methods.  Mean  values  of  p,  E* 
and  of  the  plateau  stress  a*  are  listed  in  Table  3. 

2. 2  Woods 

We  also  tested  a  number  of  woods,  chosen  to  give  a  range  of 
relative  densities  between  0.05  and  0.5.  Samples  of  well-seasoned  woods, 
roughly  20  mm  x  20  mm  x  40  mm,  were  cut  with  the  long  direction  parallel 
to  a  radius  of  the  trunk  (Fig.  7)  and  parallel  to  the  axis  of  the  trunk 
(Fig.  8).  The  samples  were  stored  for  10  days  to  reach  an  equilibrium 
moisture  content  (roughly  12  %)  and  tested  in  compression  at  18  °C  and 
a  strain  rate  close  to  10“3/s. 

3.  MODELS  FOR  THE  MECHANICAL  PROPERTIES 

When  a  cellular  material  is  compressed,  the  cell  walls  deform.  The 
deformation  modes  (bending,  buckling,  plastic  collapse  and  fracture)  are 
known  from  studies  of  model  cells  (2).  Two-dimensional  models  (shown  in 
insets  of  subsequent  figures)  can  be  analysed  accurately  (2).  With  this 
knowledge,  a  kind  of  dimensional  analysis  of  3- dimensional  cellular 
solids  (1,  3,  4)  becomes  possible.  This  analysis. 


giving; 


(6) 
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where  C2  is  a  constant.  The  shear  modulus  scales  in  a  similar 
way,  because  shear  deformation  in  a  foam  also  causes  simple  bend¬ 
ing  of  the  cell  walls.  (For  closed  cells,  I«  £t3  and  p/p  *  t/l 
giving  instead  E*/Es  *  (p/pg) 3)  • 

Data  are  compared  with  eqn.  (6)  in  Fig.  9.  The  full  line  is 
a  plot  of  eqn.  (6)  with  C2  =  1.  It  gives  a  good  description  of 
a  wide  range  of  materials  and  densities.  (We  find  that  plastic 
foams  deviate  systematically  towards  the  line  E*/E  =  (p/p  )3^2, 
because,  we  believe,  some  limited  plasticity  occurs  even  under 
small  loads.  Section  3.3,  below,  explains  the  power  of  3/2.) 

Poissons  ratios  v  have  been  measured  for  cellular  solids  (1,3).  In  the 
linear-elastic  regime, v  %  ,  although  in  the  plateau  regime  it  is  almost  zero 


3. 2  Elastic  Buckling 

Flexible  foams  show  extensive  non-linear  elasticity.  I"  is 
caused  by  the  elastic  buckling  of  the  cell  walls  (8  ,18),  as 
shown  in  the  inset  of  Fig.  10,  and  it  is  this  that  gives  the  plateau 
of  the  stress-strain  curve  for  elastomeric  foams. 


The  critical  load  at  which  a  column  of  length  l.  Young's 
modulus  Eg  and  second  moment  of  area  I  buckles,  is  given 
by  Euler's  formula: 


F 

cr 


oe 


E  I 
s 

lz 


(7) 


If  this  load  is  reached  for  a  layer  of  cells  spanning  the  section, 

they  buckle,  initiating  the  elastic  collapse  of  the  foam.  For  the 

* 

3-dimensional  open-cell  foam  of  Fig.  1  the  stress  a  . 


at  which 
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this  occurs  is  proportional  to  F^/t2.  Using  the  facts  that 

I  «  and  p/p  *  (t/l)2  we  obtain  the  elastic  collapse,  or 

s 

plateau,  stress : 


★ 


It  is  valid  for  relative  densities  below  0.3  .  At  higher  densities,  the 
cell  walls  are  too  short  and  stocky  to  buckle;  instead  theyyield  or  crush. 

it 

Data  for  for  elastomeric  foams  are  compared  with 

eqn.  (8)  in  Fig.  10.  They  are  well  fitted  by  the  equation  with 
C3  =  0.05. 

3. 3  Flastic  Collapse 

If  the  cell-wall  material  yields  plastically,  as  do  metals  and 
many  polymers,  then  the  foam  as  a  whole  shows  a  plateau  caused 
by  plastic  collapse.  It  occurs  when  the  moment  on  the  inclined 
cell  walls  exceeds  the  fully  plastic  moment,  creating  plastic 
hinges  (2,3,19,20)  as  shown  in  the  inset  of  Fig.  11.  For  a 
beam  of  square  section  of  side  t,  the  fully  plastic  moment  is: 


M  =  o  t3 

p  y 


(9) 


The  moment  is  proportional  to  FI,  and  (as  before)  the  force  F 

is  proportional  to  ol2.  Combining  these  results  with  eqn.  (5) 

* 

we  find  the  plastic  collapse,  or  plateau  stress  Opj  to  be: 


•  '4 


t\ '.! 


c4 


(10) 


Data  for  the  plateau-stress  of  plastic  foams  are  plotted  in 
Fig.  11.  They  are  well  fitted  by  eqn.  (10)  with  C4  =  0.3  for  rela¬ 
tive  densities  less  than  0.6  (1);  at  higher  densities  the  cell  edges 
are  too  short  and  stocky  to  bend  plastically  instead,  they  shear. 

3. 4  Brittle  Crushing 

Brittle  foams  (ceramics,  and  certain  rigid  polymers)  collapse 
by  yet  other  mechanisms:  brittle  crushing  in  compression  (23), 
brittle  fracture  in  tension  (24,25).  Let  the  modulus  of  rupture 
(the  maximum  surface  stress  at  the  instant  of  fracture)  for  the 
cell-wall  material  be  c^.  Then  a  cell  wall  will  fail  as  shown 
in  the  inset  to  Fig.  12  when  the  moment  acting  on  it  exceeds: 


t3 


The  moment  due  to  F  is  proportional  to  Ft,  and  the  stress  to 

F/2.2.  Combining  these  with  eqn.  (5)  gives  -the  crushing 

* 

strength  of  the  foam: 


* 


(11) 


The  limited  experimental  data,  shown  in  Fig.  12,  are  consistent  with 


Cy  =  0.65  but  are  insufficient  to  give  much  confidence  in  eqn.  (11).  But 
observations,  reviewed  elsewhere  (4),  suggest  that  the  model  has  the 
correct  physical  basis,  and  we  shall  employ  it  in  subsequent  sections. 


3. 5  Densification  and  the  Shape  of  the  Stress-Strain  Curve 

The  plateau  ends  when  the  folding  cell  walls  begin  to  touch.  During 

elastic  buckling  or  plastic  collapse,  the  foam  compresses  axially  with 

almost  no  lateral  spreading  (v  ^  0) .  Then  simple  geometry  gives  the 

relative  density  p  (e)/p  after  a  nominal  compressive  strain  e  as: 

s 


P  (c)  _  P_  f _ 1 

P  p,  1  - 
s 


■) 


e 


(12) 


where  p/pg  is  the  initial  relative  density.  Densification  is  complete, 

and  the  stress-strain  curve  becomes  almost  vertical,  when  p  (c)/p  .  =  1» 

s 

when  the  strain  is: 

Cf  =  1  -  p/ps  (13) 


We  find  experimentally  that  the  end  of  the  plateau  corresponds  to  .. 
P  (0/ps  ^  0.33  (the  solid  occupies  of  the  total  volume) when  the 
strain  is: 


es  -  1  -  3  (p/ps)  (14) 

Consider  now  the  shape  of  the  stress-strain  curve  for  flexible 
foams.  In  any  sample  of  the  foam,  there  is  a  distribution  of  cell 
edge-lengths  and  angles.  The  plateau  starts  when  a  layer  of  cells 
(those  with  longest,  or  most  favourably  oriented  edges)  buckle  (eqn.  7). 

A  small  increase  in  stress  is  needed  to  cause  more  cell  edges  (those  which 
are  slightly  shorter)  to  buckle.  We  postulate  that  the  length  of 


the  edges  which  are  about  to  buckle  after  a  strain  e,  i  (c),  is  given  by: 


rvr ■  i!  ■  i 1 


my  .v  ■  j.  .  i-  ■  .V  --  • 
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Me)  =  * 


(1  -  (£-^-)^) 
_ ps 

0  1  -  Cp/ps)^ 


C15) 


§ 


s rta 


(The  value  l  (e)  ranges  from  lQ  at  the  start  of  deformation 
to  zero  when  full  density  is  reached.)  Using  eqn.  (12)  we  obtain: 


L 


Ce)  -  Zq 


(1 


1  -  (p/ps) 


¥ 


(16) 
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The  argument  of  Section  3.2  can  now.be  repeated.  The  force  F 
on  the  cell  wall  is  related  to  iQ  throughout  the  test  by 
F  =  a  £*.  Collapse  of  cells  with  edge  length  £  (given  by 
eqn.  (16))  occurs  when  this  force  exceeds  the  buckling  load 
given  by  eqn.  (7).  Combining  these  equations  gives: 


f-  -  0.05  (2-)  < 
s  s 


1  -  (P/Ps) 


)  2 


1  -  (P/P 


s  (i  -  ty 


(17) 


which  reduces  to  our  original  eqn.  (8)  at  small  strains,  but  becomes 
infinite  (approximating  Eg)  at  the  strain  given  by  eqn.  (13) . 

A  similar  analysis  can  be  made  for  plastic  foams.  We  postulate  that 
the  length  of  the  beams  which  are  about  to  bend  plastically  decreases  as 
strain  proceeds  (because  the  long  ones  bend  first).  If  the  length  at 
a  strain  e  is  given  (as  before)  by  eqn.  (16),  then  by  the  argument 
of  Section  3.3,  the  strength  of  a  plastic  foam  follows: 


(18) 


o_ 

E. 


o  % 

«•*  s*  tf-5  < 


1  -  (p/p,)- 
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Cp/ps  (r-  e 


)?' 


This  reduces  to  eqn.  (10)  at  small  strains,  but  becomes  infinite 
(approximating  Eg)  as  the  foam  is  compressed  to  the  solid 
density  (eqn.  13). 

These  two  results  give  an  approximate  description  of  the 
stress-strain  curves  for  flexible  and  for  rigid-plastic 
foams,  in  the  fields  of  plastic  collapse  and  densification. 

They  are  combined  with  the  equations  for  linear-elasticity, 
in  the  next  section,  to  construct  deformation-mode  maps 

tovtnjvj  -  o-tsor^tCo»\  oUA^nawvs . 

4.  CONSTRUCTION  OF  DEFORMATION-MODE  MAPS 

We  have  seen  that  when  an  elastomeric  foam  is  compressed,  it  first 
deforms  in  a  linear-elastic  way;  then  its  cells  buckle  to  give  non¬ 
linear  elasticity;  and,  finally^  the  cells  collapse  completely  and  the 
stress  rises  rapidly  as  their  faces  and  edges  are  forced  together.  A 
plastic  foam  behaves  in  a  somewhat  similar  way ,  except  that  linear 
elasticity  is  now  followed  by  plastic  collapse,  and  the  ultimate 
forcing  together  of  the  cell  walls.  With  brittle  foams,  progressive 
crushing  can  again  lead  to  a  plateau  which  ends  when  the  material  is 
completely  crushed.  We  have  seen,  too,  that  each  of  these  processes 
can  be  modelled  adequately  by  using  classical  beam  theory  to  analyse 
the  deformation  of  cell  walls.  The  analysis  relates  foam  properties 
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to  the  relative  density  (p/pg)  and  to  the  properties  of  the  material 
of  which  the  foam  is  made.  The  relations,  summarised  in  Table  4, 
suggest  a  normalisation  which  brings  the  properties  of  foams  with 
the  same  relative  densities  into  coincidence.  Then  the  properties 
of  an  entire  family  of  foams  can  be  shown  as  a  deformation  mode  map 
(4),  of  which  Figs.  13  to  16  are  examples.  The  map  has  axes  of 
normalised  compressive  stress  cr/Es,  and  the  compressive  strain 
e.  It  shows  the  fields  in  which  each  mode  of  deformation  (linear 
elasticity,  non-linear  elasticity,  plastic  collapse  and  so  forth) 
is  dominant.  Superimposed  on  the  fields  are  stress-strain  con¬ 
tours  for  constant  (initial)  relative  density. 

4. 1  Elastomeric  Foams 

Figs.  13  and  14  show  mechanism-mode. maps  for  elastomeric  foams. 

Fig.  13  shows  the  experimental  stress-strain  curves  for  polymeric 
foams;  Fig.  14  is  based  on  the  theory  alone.  Mechanism  field  boun¬ 
daries  (heavy  lines)  are  shown  on  both  figures.  They  were  constructed 
as  follows. 

The  linear-elastic  regime  ends  when  elastic  buckling  begins. 

Using  eqns.  (6)  and  (8),  and  the  fact  that  c  =  Ee  in  the  linear- 
elastic  region ,  we  obtain  the  strain  corresponding  to  the  boundary  of 
the  linear-elastic  field; 

c  =  C3  =  0.05 

At  relative  densities  above  about  0.3  the  cell  walls  become  so  stocky  that 
they  can  no  longer  buckle  elastically.  The  field  boundary  thus  bends  until 
it  is  tangent  to  the  linear-elastic  loading  line  for  p/pg  =  0.3. 

Once  elastic  buckling  starts,  the  stress  is  related  to  the  strain 
by  eqn.  (17).  We  define  the  transition  from  buckling  to  densification 


as  the  line  at  which  the  relative  density  has  reached  Then,  from 
eqn.  (14)  and  (17)  the  equation  of  the  transition  is: 


Tj— •  =  5.9  x  10“ 2  (1  -  e) 2  {1  -  (^-i-^)  } 
Es  3 


It  is  plotted  as  a  heavy  line,  sloping  down  from  left  to  right  on  the  figures. 

Fig.  14  shows  a  theoretical  map.  The  contours  are  stress-strain  curves 
for  foams  of  relative  density  between  0.01  and  0.4.  They  show  a  linear 
elastic  regime  (eqn.  6)  and  a  plateau  corresponding  to  elastic  buckling; 
they  start  to  bend  upwards  when  densification  starts  (eqn.  17);  and  they 
approach  a  limiting  slope  of  Eg  when  densification  is  complete  (eqn.  13). 
Within  the  field  of  elastic  buckling  the  material  exists  in  two  states  at 
almost  the  same  stress  (the  linear-elastic  state  and  the  densified  state) ; 
it  is  like  the  p-V  response  of  an  ideal  gas  (or  the  temperature-entropy 
diagram  for  steam)  in  which  gas  and  liquid  states  can  co-exist.  The 
material  deforms  by  the  formation  of  densified  bands  which  thicken ,  at 
constant  stress.,  as  the  strain  is  increased,  until  the  entire  material 
has  reached  the  dense  state. 

The  figure  describes  the  overall  response  of  all  isotropic,  flexible 
foams  in  compression.  In  tension,  flexible  foams  are  roughly  linear- 
elastic  to  rupture. 


Rigid.  Plastic  Foams 


Plastic  foams,  like  the  elastic  ones,  show  three  regions:  linear 
elasticity,  plastic  collapse  and  densification  -  though  now  the  strain 
beyond  the  linear-elastic  regime  is  not  recoverable.  Figs.  15  and  16 
are  a  pair  of  maps,  one  showing  our  experimental  stress-strain  curves  for 
plastic  foams,  the  other,  based  on  the  theory  alone.  Mechanism  field- 


boundaries  are  superimposed  on  the  stress-strain  curves. 


The  boundary  of  the  linear-elastic  field  (heavy  line)  is  obtained 
from  the  equation  for  linear  elasticity  (eqn.  6)  and  that  for  plastic 
collapse  (eqn.  10);  its  equation  is: 


^■(o4)  4 


In  constructing  the  map  we  have  taken  to  be  0.10  (a  fairly  typical 

value).  Next  to  the  linear-elastic  field  is  the  field  of  plastic  collapse. 

As  before,  two  states  of  strain  co-exist  at  almost  the  same  stress,  so  that 
complete  collapse  of  part  of  the  structure  can  occur  while  the  rest  is  still 
elastic;  the  bands  of  dense  material  broaden  with  increasing  strain.  We  de¬ 
fine  the  transition  from  collapse  to  densification  as  the  line  at  which  the 
relative  density  has  reached  Then,  from  eqn.  (14)  and  eqn.  (18),  the 
equation  of  the  transition  line  is: 


— -  =  0.19 
Es 


*  ! 
e)  (1  -  (- 


It  is  plotted  as  a  heavy  line  sloping  down  from  left  to  right  on  Figs.  15 
and  16. 

Fig.  16  shows  a  theoretical  map  for  plastic  foams  with  o y/^s  ~  0.10. 

It  shows  fields  of  elastic  deformation,  plastic  collapse  and  densification. 
Superimposed  on  the  fields  are  stress-strain  curves  for  foams  with  densities 
from  0.01  to  0.4. 

The  figure  shows  the  overall  response  of  isotropic,  plastic  foams  in 
compression.  It  is  less  general  than  the  map  for  elastomeric  foams  because  it 
must  be  constructed  for  a  particular  value  of  °y/Es*  But  the  equations  show 
that  the  boundaries  are  not  very  sensitive  to  its  value,  and,  for  a  given 
material,  the  diagram  shows  the  behaviour  for  all  densities. 


The  behaviour  of  plastic  foams  in  tension  resembles  that  in  com¬ 


pression,  truncated  by  fracture. 

4.  3  Brittle  Foams 

Rigid  foams  show  linear-elastic  behaviour  to  fracture.  In  compres¬ 
sion,  the  foam  crushes  at  constant  stress  (eqn.  11),  and  since  the  crush- 

« 

ing  equation  has  the  same  form  as  that  for  plastic  collapse,  the  behaviour 
will  resemble  that  of  Figs.  IS  and  16.  If  the  foam  is  contained,  it  will 
densify  at  the  strain  given  approximately  by  eqn.  (18),  with  a y/Eg  re¬ 
placed  by  Of! Eg. 

In  tension,  linear  elastic  behaviour  is  truncated  by  fast,  brittle  fracture 
■The  fracture  mechanics  of  foams  (28)  need  not  concern  us  here. 

4. 4  Woods 

These  ideas  can  be  applied,  in  an  approximate  way,  to  the  compressive 
deformation  of  wood.  Woods  are  cellular  solids,  composed  of  mixed  polymers 
(cellulose,  lignin,  hemicellulose) ,  and  with  a  relative  density  ranging  from 
less  than  0.0S  (balsas)  to  almost  1  (lignum  vitae).  Compressed  across  the 
grain,  wood  behaves  like  a  rigid-plastic  foam.  Stress-strain  curves  for 
the  woods  we  tested  are  shown  in  Fig.  17,  plotted  on  axes  of  cr/Es  and 
e  (data  in  Tables  2  and  3). 

The  field  boundaries  of  Fig.  17  were  constructed  as  follows.  In 
radial  compression,  Young's  modulus  for  woods  is  found  to  follow  the  same 
law  as  elastic  foams  (29)  : 

E*  =  10  (£-)2  GN/m2  (22) 

Ps 

Compressive  collapse  in  the  radial  direction  eqn.  (29)  starts  at  the  stress: 

2 

o*  *  135  (£- )  MN/m  2 

V 


(23) 


1? 


if 

V-v 

.  i . 

*  v. 

•y  e  =  0.004 


When  woods  are  compressed  axially,  they  collapse  by  a  process  which 
involves  the  kinking  (50)  or  fracture  (29,30)  of  the  cell  walls. 

This  occurs  at  essentially  constant  stress  until  a  critical  state 
is  reached,  when  the  stress-strain  curve  bends  sharply  upwards. 
Experimentally,  we  find  that  the  bend  upwards  occurs  at  a  strain  ec 
corresponding  to  a  density  of  roughly  0.5: 

■ 1  - 2 

s 

Combining  this  with  eqn.  (27)  gives  the  field  boundary: 


£=!>  10-3  (1  -  c) 

Es 

It  is  plotted  as  a  heavy  line  on  Fig.  18.  The  figure  summarises 
the  compressive  behaviour  of  woods  in  the  axial  direction.  Both  dia¬ 
grams  allow,  by  interpolation,  or  the  use  of  the  equations,  the  approxi¬ 
mate  prediction  of  the  stress-strain  curve  for  woods  of  other  densities. 


5.  CONSTRUCTION  OF  ENERGY-ABSORPTION  DIAGRAMS 

The  commonest  use  of  foams  is  in  packaging.  The  aim  is  to  absorb 
energy  (usually  the  kinetic  energy)  of  the  packaged  object  when  it  is 
dropped,  or  is  accelerated  or  decelerated  in  some  other  way,  while  at 
the' same  time  keeping  the  force  on  the  object  below  the  limit  which 
will  cause  damage.  In  selecting  a  material  for  the  package,  we  need 
to  know  the  energy  that  can  be  absorbed  without  the  stress  exceeding 
a  critical  value.  We  show  below  that  there  is  an  optimum  foam  density 
for  a  given  package.  If  the  density  is  too  low,  the  foam  'bottoms  out' 
(with  a  sharp  increase  in  stress)  before  enough  energy  has  been  absorbed.  If 


it  is  too  dense,  the  stress  exceeds  the  critical  value  before  enough 
energy  has  been  absorbed. 


Figure  19  shows,  inset,  part  of  a  stress-strain  curve  for  a  foam. 


It  is  linear-elastic  to  e  ,  and  thereafter  follows  the  stress-strain 

o 

a  (e)  curve  described  approximately  by  eqn.  (17)  or  (18).  The  area,  up 
to  the  strain  e*,  where  the  stress  is  a,  (shaded  on  Fig.  19)  is: 

.e0  e* 

a  dc  +  |  o(e)  de 
o  e0 


f  f 

W  =  a  de  + 


or 


JL  =  I_Le2  +  f 

Es  2  Es  0  j 


o(e) 


de 


This  equation  was  integrated  numerically,  using  eqn.  (17)  for  a(e)  for 
elastic  foams  and  eqn.  (18)  for  plastic  foams t  to  give  Figs.  19  and  21. 

Similar  diagrams  can,  of  course,  be  constructed  directly  from  its  experi¬ 
mental  stress-strain  curves  (such  as  those  shown  in  Figs.  3  to  6)  by 
measuring  the  area  W/Es  up  to  the  strain  e*  corresponding  to  the 
stress  o/Es.  Such  diagrams  are  shown,  for  comparison,  in  Figs.  20  and  22. 

The  energy  diagram,  calculated  for  elastic  foams  (Fig.  19)  shows  the 
normalised  energy  absorbed  per  unit  volume  of  foam,  W/Eg,  plotted  against 
the  peak  stress  o/Es,  for  a  range  of  densities  p/pg.  Normalised  in  this 
way,  the  diagram  describes  all  elastic  foams.  If  the  critical  damage 
stress  is  selected,  then  the  diagram  gives  the  foam  density  which  will  absorb 
the  greatest  amount  of  energy  without  this  stress  being  exceeded.  The  lower 
the  peak  stress,  the  lower  is  the  optimum  foam  density.  As  an  example,  using 
Fig.  19  or  20,  the  optimum  foam  density  for  a  critical  damage  stress  cr/Es  = 
is  p/pg  =  0.1..  Choosing  the  right  density  is  important:  it  can  easily  give 
a  factor  of  10  greater  absorption  of  energy  than  a  wrongly-chosen  foam. 


Figures  21  and  22  show  energy  diagrams  for  a  plastic  foam  with 
a  /E  =0.1.  Such  diagrams  aTe  less  general  than  those  for  elastic 

y  * 

foams,  because  a  / E  must  be  specified;  but  a  single  diagram  still 

y  s 

describes  all  foams  made  of  a  given  material  (polystyrene,  for 
example).  Here,  too,  there  is  an  optimum  foam  density  for  a  given 
energy  absorption  and  peak  stress.  The  diagram  allows  it  to  be 
chosen. 

The  figures  show  that,  for  a  given  material  (and  thus  Es) ,  there  is 
a  maximum  energy  which  can  be  absorbed  for  a  given  peak  stress.  The 
envelope  shown  as  a  broken  line  on  Figs.  19  and  21  divides  the  diagram 
into  an  accessible  region  (below  the  line)  and  an  inaccessible  one  (above). 
For  elastomeric  foams,  the  equation  of  the  line  is  approximately: 

7/8 

W/E  =  0.11  {§-) 

3 

It  applies  to  all  elastic  foams.  For  the  plastic  foam  with  a  / E  =0.1, 

y  s 

the  equation  of  the  line  is  approximately: 

2/3 

W/E  =0.05  (£-) 
s  Es 

The  equations  from  p/ps  =0.01  to  p/pg  =  0.3. 

6.  CONCLUSIONS  AND  APPLICATIONS 

6. 1  Cone lusions 

When  a  cellular  solid  with  a  relative  density  below  0.30  is  compressed, 
it  shows  a  stress-strain  curve  with  three  parts:  a  linear-elastic  part, 
a  long  plateau,  and  a  regime  of  final  densification.  The  properties  of 
an  entire  class  of  such  solids  can  be  summarised  as  a  deformation -mode 
map,  which  shows  how  each  of  the  three  parts  changes  as  the  relative 


r-w 


t 


density  changes.  The  method  can  be  applied  to  man-made  foams  (Figs.  13  to 
16)  and  to  natural  materials,  such  as  woods  (Figs.  17  and  18),  to  give 
diagrams  which  summarise  the  stress -strain  response  of  each  class  of 
material. 

The  mechanical  properties  of  cellular  solids  can  be  modelled  with 
precision.  The  models  lead  to  constitutive  laws  (stress-strain  relations) 
which  have  been  thoroughly  tested  in  simple  compression.  They  are  summa¬ 
rised  in  Table  4.  The  constitutive  laws  can  be  used  to  construct  deforma¬ 
tion-mode  maps,  and  can  be  integrated  to  construct  energy-absorption  dia¬ 
grams.  These  diagrams  (Figs.  19  to  22)  show  how  the  energy  absorbed, 
per  unit  volume  of  foam,  depends  on  the  density  of  the  foam  and  on  the 
stress.  The  diagrams  show  that  there  is  an  optimum  foam  density  for  a 
given  packaging  or  energy-absorbing  application. 

6.2  Design  with  Cellular  Solids 

The  equations  derived  in  the  text,  and  summarised  in  Table  4,  pro¬ 
vide  the  basic  information  for  design  with  foams  in  load-bearing  applica¬ 
tions.  The  tests  described  above  all  involved  simple  compression. 

Under  multiaxial  loads,  the  behaviour  is  more  complicated.  At  small 
strains  (e  <  5  %) ,  while  the  material  is  linear-elastic,  it  behaves 
like  any  other  elastic  solid  (1,3)  with  Young's  modulus  E*  given  by 
eqn.  (6)  and  Poisson's  ratio  v  =  -$■.  But  once  the  plateau-stress  has 
been  reached,  the  behaviour  changes:  the  extensive  deformation  at  nearly 
constant  stress  involves  a  large  volume  change,  but  almost  no  lateral 
strain  (1,3,26,27)  so  that  v  ^  0 .  Then  the  material  deforms  under  a 
multiaxial  state  of  stress  when  the  maximum  principal  stress  (not  the 
octahedral  shear  stress)  reaches  the  critical  value  a*  (eqn.  8  for  flex¬ 
ible  foams,  eqn.  10  for  plastic  foams).  Because  of  this,  the  indentation 


as  in  dense 


hardness  of  plastic  foams  is  equal  to  o*^  (not  3o*^ 
solids)  and  the  force  needed  to  compress  a  foam  is  the  same  whether 
it  is  free  at  its  edges  or  constrained  there  (15)  .  Once  densification 
starts,  of  course,  the  properties  revert  towards  those  of  conventional 
solids. 


The  deformation  mode  map  for  a  class  of  foams  provides  a  compact 
summary  of  the  mechanical  response  of  the  class,  in  particular,  the 
maps  give  a  rational  way  of  selecting  the  material  and  density  which 
will  give  a  desired  stress-strain  response.  Similarly,  the  energy- 
absorption  diagram  for  a  class  of  foam  summarises  the  energy  absorbing 
capacity  of  all  members  of  the  class.  The  diagrams  identify  the  optimum 
foam  material  and  density  for  a  given  packaging,  padding  or  cushioning 
application. 
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TABLE 
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EXPERIMENTA 


SAMPLE 


p  (kt|/m3) 


El ASTOMFRIC  FOAMS 


Dunlop 


F relen 


D1 
D14 
D17 
F  30 
F70 
F 1 20 
FI  75 


1i|.4 

32.4 
51.7 

29.4 
69.2 

120.0 

138 


F  2  50  \ 

360 

PLASTIC  FOAMS 

Rotucel  1  31 

34 

51 

51 

.6 

71 

70.4 

no 

124 

170 

160 

190 

165 

.7 

brittle  foams 

SI 

100 

S2 

200 

S3 

280 

S4 

320 

WOODS 

Radial 

— 

Axial 

Balsa 

103 

103 

Balsa 

124 

124 

Willow 

403 

355 

Pine 

443 

- 

Willow 

S20 

385 

Mahogany 

615 

- 

Beech 

111 

73S 

Beech 

737 

T  ’  *  '  ~ 

764 

1 1  *  *i 

MEASUREMENTS 


£* (MPa ) 

<s  (MPa 

C.  054 

0.002 

0.062 

0.0025 

0.057 

0.005 

0.275 

0.04 

1.10 

0.1 

3.82 

0.24 

4.58 

0.3 

23.46 

2.56 

Tension 


43 

73.4 

96.2 

194 

280 

432 


Compr . 


20.8 

26.1 

56.9 

129.7 

196.7 
258.9 


0.4 

0.8 

1.6 

3.0 

5.2 

6.0 


8-50 

6-22 

2-8 

25-44 


0.02  -0.18 
0.035-0.16 
0.072-0.21 
0.093-0.16 


Radial  1 

Axial 

1  1  — 

Radial 

Axial 

31 

42D 

0.80 

8.0 

79 

258 

1.20 

7.0 

198 

1270 

3.70 

22.0 

647 

- 

6.50 

- 

321 

2203 

8.2S 

22.5 

610 

“  1 

- 

- 

1160 

7393 

• 

50.0 

1569 

7500 

- 

50.0 

TABLE  4;  EQUATION  FOR  STIFFNESS  AND  STRENGTH  OF  CELLULAR  SOLIDS 


Property 


Equation 


Young's  modulus 


—  =  (*-) 2 
Es  V 


Plateau  stress, 
flexible  foams 


=  0.05  <S_)1 
0  - 


Plateau  stress, 
plastic  foams 


a  i 

-Ei  =  0.3  (E-) 
0„  Pe 


Approximate  plateau  -  0  65  fp  V 

stress,  brittle  foams  a£  Lp„J 


r  1  -  (£-)*  32 

Stress-strain  response,  o_  n  ,p  -,2 _ ps  _ 

flexible  foam  Eg  ~  ‘  lpg'  j  _  ^p__  ^ _ ——jfc 

ps  1  *  e  ■ 


Stress-strain  response,  | 
plastic  foam 


■f-  -  °.3  (M 
°y  ps 


*  1  -  (r> 1  | 

pp)  l1  •  t  <rh^j 
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FIGURE  CAPTIONS 


Fig.  1.  Schematic  of  cellular  material,  showing  dimensions. 


Fig.  2.  Microstructures  of  the  four  types  of  foam: 

(a)  polyethylene  p/pg  =0.115;  Cb)  polyurethane  p/p  =0.029 

(c)  polymethylacrylimide  p/p  =0.105;  (d)  mullite  p/ps  =0.062 

s  s 

Fig.  3.  Compressive  stress- strain  curves  for  flexible  polyethylene  foams. 


Fig.  4.  Compressive  stress-strain  curves  for  flexible  polyurethane  foams. 


Fig.  5.  Compressive  stress-strain  curves  for  rigid  (plastic)  polymethyl- 
acrylimid  foams. 


Fig.  6.  Compressive  stress-strain  curve  for  a  brittle  mullite  foam. 


Fig.  7.  Compressive  stress-strain  curves  for  woods,  in  the  radial 
direction  (perpendicular  to  the  grain) . 


Fig.  8.  Compressive  stress-strain  curves  for  woods  in  the  axial 
direction  (parallel  to  the  grain) . 


Fig.  9.  The  relative  Young's  modulus,  E*/Es,  plotted  against  relative 
density.  p/ps. 


Fig.  10.  The  relative  elastic  collapse  stress,  a|l/Es,  plotted 
against  relative  density,  p/ps. 


Fig.  11.  The  relative  plastic  collapse  stress,  Cpi/Oy,  plotted 
against  relative  density  p/ps.  '  3 


Fig.  12.  The  relative  crushing  strength.  Of/of,  plotted  against 
relative  density  p/ps,  for  brittle  foams. 


Fig.  13.  A  deformation  mode  map  for  flexible  foams.  It  shows  the 
data  of  Figs.  4  and  5,  normalised.  The  construction  of 
the  field  boundaries  is  described  in  the  text. 


Fig.  14.  A  deformation-mode  map  for  flexible  foams,  constructed 
entirely  from  the  equations  developed  in  the  text. 


Fig.  15.  A  deformation-mode  map  for  plastic  foams.  It  shows 

the  data  of  Fig.  6,  normalised.  The  construction  of  the 
field  boundaries  is  described  in  the  text. 


Fig.  16.  A  deformation-mode  map  for  plastic  foams,  constructed 
entirely  from  the  equations  developed  in  the  text. 

Fig.  17.  A  deformation-mode  map  for  woods,  tested  in  the  radial 
direction  (across  the  grain) . 

Fig.  18.  A  deformation-mode  map  for  woods  tested  in  the  axial 
direction  (along  the  grain) . 


Fig.  19.  An  energy- absorption  diagram  for  elastomeric  foams, 
constructed  from  the  equations  given  in  the  text. 

The  broken  line  divides  the  diagram  into  an 
accessible  and  an  inaccessible  region. 

Fig.  20.  An  energy- absorption  diagram  for  elastomeric  foams, 
constructed  by  measuring  the  areas  under  the  stress- 
strain  curves  of  Figs.  3  and  4.  It  is  directly  com¬ 
parable  with  Fig.  19. 


Fig.  21.  An  energy -absorption  diagram  for  plastic  foams,  const¬ 
ructed  from  the  equations  given  in  the  text,  with 
Oy/Es  =  0.1.  The  broken  line  divides  the  diagram 
into  an  accessible  and  an  inaccessible  region. 

Fig.  22.  An  eneTgy-absorption  diagram  for  polymethacrylimid 
foams,  constructed  by  measuring  the  areas  under  the 
stress-strain  curves  of  Fig.  5.  It  is  directly  com¬ 
parable  with  Fig.  21. 
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Schematic  of  cellular  material,  showing  dimensions. 
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Fig.  2.  Microstructures  of  the  four  types  of  foam: 

(a)  polyethylene  p/p  =  0 . 1 1  5 ;  (b)  polyurethane  p/p  =0.029 

(c)  polymcthylacrylimidc  p/p*  =0.103;  (d)  mullite  p/pS  =0.062 
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Compressive  stress- strain  curves  for  flexible  polyurethane  foams 


Fig.  7.  Compressive  stress-strain  curves  for  woods,  in  the  radial 
direction  (perpendicular  to  the  grain) . 


Fig.  8.  Compressive  stress-strain  curves  for  woods  in  the  axial 
direction  (parallel  to  the  grain). 
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Fig.  10.  The  relative  elastic  collapse  stress,  o<» 
against  relative  density,  p/ps. 
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The  relative  plastic  collapse  stress,  cr*i/c^,  plotted 
against  relative  density  p/ps* 
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CRUSHING  STRENGTH /MODULUS  OF  RUPTURE 


The  relative  crushing  strength,  of/of,  plotted  against 
relative  density  p/ps,  for  brittle  foams. 


COMPRESSIVE  STRAIN  € 


Fig.  13.  A  deformation  mode  map  for  flexible  foams.  It  shows  the 
data  of  Figs.  4  and  5,  normalised.  The  construction  of 
the  field  boundaries  is  described  in  the  text. 
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Fig.  14.  A  dcformation-modc  map  for  flexible  foams,  constructed 
entirely  from  the  equations  developed  in  the  text. 


Fig.  15.  A  deformation-mode  map  for  plastic  foams.  It  shows 

the  data  of  Fig.  6,  normalised.  The  construction  of  the 
field  boundaries  is  described  in  the  text. 
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COMPRESSIVE  STRAIN,  € 

Fig.  16.  A  deformation-mode  map  for  plastic  foams,  constructed 
entirely  from  the  equations  developed  in  the  text. 


Fig.  19.  An  energy- absorption  diagram  for  elastomeric  foams, 
constructed  from  the  equations  given  in  the  text. 
The  broken  line  divides  the  diagram  into  an 
accessible  and  an  inaccessible  region. 


An  energy-absorption  diagram  for  elastomeric  foams, 
constructed  by  measuring  the  areas  under  the  stress- 
strain  curves  of  Figs.  3  and  4.  It  is  directly  com¬ 
parable  with  Fig.  19. 
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ABSTRACT 

This  paper  explores  the  possibility  of  constructing 
deformation-mechanism  maps  for  amorphous  polymers.  Five  regimes  are 
identified:  the  glassy  regime,  the  visco-elastic  regime,  the  rubbery  regime 
and  the  regime  of  viscous  flow  (melting),  truncated  by  decomposition. 
Constitutive  equations  for  each  regime  are  assembled  and  adapted  ;o  give  a  good 
description  of  a  large  body  of  experimental  data  for  amorphous 
polymethylmethacrylate  and  polystyrene.  The  adjusted  laws  are  then  used  to 
construct  diagrams  which  relate  the  time- and  temperature-dependent  modulus, 

E  (t,T),  to  the  temperature  and  the  loading  time  (or  frequency'.  The 
diagrams  are  divided  into  fields  corresponding  to  the  five  regimes. 

A  diagram  summarises  the  small-strain  mechanical  behaviour  of  the  polymer 
over  a  wide  range  of  conditions. 
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TABLE  1:  SYMBOLS,  DEFINITIONS  AND  UNITS 

Temperature  (X) 

Tensile  or  compressive  stress  (MPa) 

Tensile  or  compressive  strain 

Tensile  or  compressive  strain  rate  (s-1) 

Time  ( s ) 

.  — i  . 

Frequency  (s  ) 

Time  and  temperature  dependent  Young's  modulus  (MPa) 

Modulus  at  0  K  (MPa) 

Temperature  coefficient  of  modulus 
Glass  temperature  (K) 

Depolymerisation  temperature 
Molecular  weights  (kg/ mol) 

Mean  activation  energy  (kJ/mol) 

Fractional  standard  deviation  of  activation  energy 
Fractional  modulus  drop  at  a  transition 
Shift  factor  for  time-temperature-equivalence 
Viscosity  (Ns/m* ) 

Pre-exponential  for  viscosity  (Ns/m*  ) 

Bulk  thermal  expansion  coefficient  below  and  above  T^  (K~l) 
W.L.F.  constants  (see  eqn.  15) 

Density  (kg/m3  ) 

Fractional  free  volume 

Internal  energy  per  unit  volume  (J/m3  ) 

Gas  constant  (3.314  J/mol  X) 

Energy  of  formation  of  vacancies  in  a  crystal  (J/mol) 

Energy  of  motion  of  vacancies  in  a  crystal  (J/mol) 

A  pre-exponential  frequency  factor  (s~l) 


1. 


INTRODUCTION 
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When  a  polymer  is  loaded,  it  suffers  deformation  which,  in  general, 
increases  with  time  of  loading.  For  uniaxial  loading,  the  resistance  to 
small-strain  deformation  is  conventionally  measured  by  the  time-and 
temperature- dependent  modulus,  E  (t,T)  (from  now  on  simply  called  E). 
If,  for  instance,  a  constant  stress  o  is  applied  to  a  sample  of  the 
polymer,  giving  a  strain  e  (t,T)  after  a  time  t  at  temperature  T, 
then: 


Linear  amorphous  polymers  like  polymethylmethacrylate  ( PMMA)  or 

polystyrene  (PS)  show  five  distinct  regimes  of  deformation,  in  each  of  which 

the  modulus  has  certain  characteristics,  illustrated  in  Fig.  1. 

1  a.)  The  glassy  regime,  with  a  modulus  of  between  1  and  10  GPa,  associated 
with  stretching  and  bending  of  intermolecular  bonds,  and  showing 
only  a  slight  time  dependence  associated  with  a  number  of 
secondary  relaxations . 

(b)  The  glass-transition  regime,  in  which  the  modulus  drops  steeply  from 
around  1  GPa  to  near  1  MPa  with  increasing  temperature  or  time  of 
loading . 

(c)  The  rubbery  regime,  with  a  modulus  of  between  1  and  10  MPa,  associated 
with  the  rubber-like  sliding  of  the  long-chain  network  of  molecules, 
constrained  by  entanglements  which  behave  (physically)  like 
cross-links . 

(d)  The  viscous  regime,  at  temperatures  well  above  the  glass  transition 
temperature,  in  which  the  polymer  can  be  thought  of  as  a  viscous 
liquid;  its  molecules  move  relative  to  each  in  a  snake-like  manner 
(reptation)  which,  when  biased  by  stress,  leads  to  viscous  flow. 

*'e)  A  regime  of  decomposition  in  which  chemical  breakdown  begins. 

Each  regime  is  associated  with  a  certain  range  of  modulus.  Each  has 
been  studied  and  modelled  in  more  or  Less  detail,  and  constitutive  laws 
for  each  have  been  formulated  -  they  are  reviewed  below.  Experiment-’ 1 
data  for  E,  in  each,  are  available  for  PMMA  and  PS.  Our  aim 


here  is  to  fit  the  constitutive  laws  to  the  data  (requiring  some 
modification  to  the  models)  and  to  examine  how  the  laws  can  be  combined  to 
give  a  self-consistent  description  of  small-strain  deformation  from  0  °K  to 
above  the  melting  temperature,  and  for  all  practical  loading  times.  The 
results  are  assembled  into  diagrams  which  summarize  the  small-strain 
deformation  behaviour  of  each  polymer  as  a  function  of  temperature ,  time  and 
frequency . 

We  now  examine  the  five  regimes  in  more  detail,  reviewing,  and 
selecting,  among  the  constitutive  laws  for  each.  Symbols  and  their 
definitions  are  listed  Table  1. 

2.  DEFORMATION  MECHANISMS 

The  Glassy  Regime  and  the  8  ,  y  and  5  Relaxations 

Well  below  the  glass  temperature,  T  ,  linear  amorphous  polymers  have 
Young's  moduli  of,  very  roughly,  3  GPa.  This  is  a  direct  reflection  of  the 
low-stiffness  Van  der  Waals  bonds  which  bind  one  chain  to  another  as  shown 
as  dotted  lines  in  Fig.  2:  when  the  polymer  is  deformed,  it  is  these  bonds 
which  stretch  and  bend.  The  covalent  C-C  bonds  which  form  the  chain 
backbone  (full  lines)  are  about  100  times  stiff er  than  the  Van  der  Waals 
bonds,  and  their  stretching  and  bending  contributes  nothing  significant  to 
the  elastic  deformation.  Rotation  about  a  C-C  bond  is  another  matter:  the 
single  C-C  bond  rotates  so  easily  that  it  is  problems  of  steric  hindrance, 
and  thus  the  Van  der  Waals  bonds,  again,  which  limit  its 
extent  (Bowden,  1968;  Yannas ,  1974;  Yannas  and  Luise,  1982). 

If  the  internal  energy  of  the  polymer,  per  unit  volume,  is  V  (ej, 
then  the  glassy  modulus  at  0  K  is  calculated  (in  principle)  by  forming  the 
second  derivative  of  V  (e)  with  respect  to  e: 

d1  V  (e) 

- - - -  (2). 

The  difficulty  is  that  of  modelling  V  (e)  with  enough  precision  to  give 
more  than  an  order-of -magnitude  estimate  of  the  modulus.  To  do  so  requires 
a  potential  function  describing  the  Van  der  Waals  bonds.  Part  of  the 
elastic  strain  comes  from  bond  stretching;  then,  commonly,  a  Lennard—J ones 
potential  is  summed  over  all  molecules  to  give  V  (e).  Strain  can  also 


result  from  rotation  of  segments  of  a  single  molecule  about  a  C-C  bond 
such  a  way  that  the  molecule  changes  shape,  (Bowden,  1968;  Yannis,  1974). 
Then  a  potential  which  includes  steric  hindrance  (Flory,  1969)  must  be  used 
to  calculate  V  (e).  But  the  fundamental  understanding  of  intermolecular 
bonding  in  polymers  is  still  too  poor  to  allow  V  (e)  to  be  computed 
accurately  from  first  principles  in  either  case  (in  practice,  experimental 
data  for  E  are  used  to  calibrate  the  potential  functions ,  rather  than  the 
other  way  round).  Instead,  the  modulus  at  absolute  zero,  E  ,  is  obtained 
by  extrapolating  experimental  measurements  of  E  at  slightly  higher 
temperatures  to  absolute  zero.  It,  and  other  material  properties,  are 
listed  in  Table  2. 


Increasing  the  temperature  has  two  distinct  effects.  First,  thermal 
expansion  increases  the  molecular  separation  and  lowers  the  Van  der  Waals 
restoring  forces:  this  gives  a  slow  drop  in  modulus,  but  does  not  introduce 
a  time  or  frequency-dependence.  Second,  the  thermal  energy  of  the  molecules 
permits  thermally-activated  local  rearrangements  (usually,  rotations  about  a 
C-C  bond),  giving  time-dependent  strain  and  an  associated  drop  in  modulus. 

The  first  effect  leads  to  a  roughly  linear  decrease  in  modulus  with 
increasing  temperature.  Yannas  and  Luise  (1982)  develop  this  idea:  thermal 
expansion  increases  the  inter-chain  distance,  reducing  the  force  required 
for  intermolecular  deformation  by  bond  stretching.  The  final  form  of  their 
equation  is  cumbersome  but  simplifies  to  a  linear  relationship  at  low 
temperatures.  Van  Krevelen  (1976)  prefers  the  empirical  relation: 


E  =  E  /(I  ♦  2  T/T  ) 
0  g 


which,  of  course,  is  linear  at  low  temperatures,  but  drops  faster  near  T 
In  this  paper  we  use  a  linear  temperature  dependence: 


E 


(1  - 


T  * 

“m  T~’ 
m  T 

g 


(3) 


to  describe  the  drop  in  modulus  caused  by  thermal  expansion.  The 
dimensionless  coefficients  which  best  described  the  low  temperature  data 
for  each  polymer  are  listed  in  Table  2. 


The  second  contribution  comes  from  the  small  thermally-activated 


rearrangements  of  side  groups  or  short  segments  of  chain  at  "loose  sites"  in 
the  amorphous  structure  (Fig.  3).  These  relaxation  processes,  all  occurring 

at  temperatures  below  T  (and  known,  in  order  of  decreasing  temperature, 

£ 

as  the  3,  Y  and  5  relaxations)  give  additional  strain,  and  they  thereby 
reduce  the  modulus  a  little.  Each  relaxation  is  positioned  for  a  given 
loading  time  or  frequency,  about  a  characteristic  temperature  at  which 
thermal  energy  becomes  sufficient  to  activate  the  rearrangement. 

Like  similar  relaxations  in  crystalline  solids,  the  response  of  the 
material  to  load  is  conveniently  described  by  a  rheological  model.  The 
simplest  that  is  of  any  real  use  is  the  standard  linear  solid  (Fig.  4a) , 
composed  of  two  springs  of  modulus  E  and  AE,  and  a  dashpot  of  viscosity: 

n»noexp£T  (4) 

where  Q  is  the  activation  energy  for  the  process ,  R  the  gas  constant  and 

n  a  constant  which  is  chosen  to  fix  the  position  of  the  drop  in  modulus, 
o 

The  standard-linear  model  gives  a  good  phenomenological  description  of  most 
anelastic  relaxations  in  crystals,  and  relaxations  in  glasses  do,  very 
roughly,  have  the  same  characteristics.  But  more  careful  examination  of 
data  for  amorphous  polymers  shows  that  the  standard  linear  solid  is  too 
simple.  Its  response  leads  to  an  almost  step  like  drop  in  modulus  at  the 
characteristic  temperature;  the  real  relaxation  is  broader  (Fig.  5).  This 
is  no  surprise.  In  PMMA ,  for  instance,  the  S-relaxation  is  thought  to  be 
caused  by  a  motion  of  the  ester  side-group,  the  Y-relaxation  by  the  motion 
of  one  of  the  two  methyl  groups  and  the  ^-relaxation  by  the  motion  of  the 
other  one  (Fig.  3) .  The  amorphous  chain-packing  grips  some  of  these  more 
tightly  than  others,  so  that  each  relaxation  has  a  spectrum  of  activation 
energies.  The  response  is  then  more  realistically  described  by  the  parallel 
coupling  of  springs-and-dashpot  units  shown  in  Fig.  4(b),  each  dashpot 
describing  a  part  of  the  activation  energy  spectrum. 

This  arrangement  of  units,  rather  than  a  more  general  one  involving  units 
in  series  and  parallel,  is  justified  in  the  following  way.  Relaxation  involves 
the  motion  of  isolated  side  groups  -  those  which,  at  a  given  instant,  are  less 
tightly  gripped  than  the  rest.  Each  can  be  thought  of  as  a  small  spherical  volume 
in  which  viscous  deformation  takes  place,  embedded  in  an  elastic  matrix.  The 
constitutive  equation  for  a  material  consisting  of  viscous  inclusions  in  an 
elastic  matrix  is  discussed  by  3rown  (1982):  provided  the  dispersion  of 
viscous  spheres  is  dilute  (meaning  that  they  are  well  separated)  ,  the 


material  behaves  like  a  standard  linear  solid  [as  in  Fig.  4a).  If  each  of  the 
spherical  volumes  has  a  different  activation  energy  (as  is  the  case  for  the 
side-group  relaxations  discussed  here)  the  same  reasoning  leads  to  the 
conclusions  that  the  material  behaves  in  a  way  described  by  Maxwell 
elements  in  parallel,  as  in  Fig.  4b. 

As  explained,  the  spectrum  of  activation  energies  arises  from 
differences  in  local  packing  density  which  makes  it  harder  to  move 
some  side-groups  than  others.  We  have  used  a  Gaussian  distribution  of 
activation  energies  with  a  mean  Q  and  a  standard  deviation  AQ.  The 
proportion  of  units  (that  is,  the  fraction  of  the  spheres)  which  relax 
with  activation  energies  between  Q  and  dQ  is  f  (Q)  dQ  where: 


0-0. 


f  (Q)  =  -  exp  -  Vz  (- 

/2  < 


This  defines  a  weighting  function  for  each  unit  (Fig.  6)  so  that  the  modulus 
drop  associated  with  full  relaxation  of  the  i  unit  is: 


5E.  *  4E  f  (0)  dQ 


where  AE  is  the  total  modulus  drop  associated  with  the  relaxation.  This 

method  introduces  only  one  new  variable,  the  standard  deviation,  AQ,  into 

the  calculation  of  the  relaxation;  it  is  chosen  so  that  the  width  of  the 

relaxation  (Fig.  5)  matches  experiment.  Data  for  AE,  Q  and  aQ /Q 

m  m 

PMMA  and  PS  are  listed  in  Table  2. 

The  Glass-Rubber  Transition  (The  Visco-Elastic  Regime) 


As  the  temperature  is  raised,  the  Van  der  Waals  bonds  start  to  melt. 
Then  segments  of  the  previously  elastically-bonded  chains  undergo  larger 
sliding  movements  relative  to  each  other  and  the  material  behaves  in  a 
visco-elastic  way.  Within  this  regime  it  is  found  that  the  modulus  E  at 
one  temperature  can  be  related  to  that  at  another  by  a  change  in  the  time 
scale  only,  that  is,  there  is  an  equivalence  between  time  and  temperature. 


this  means  that  the  curve  describing  the  modulus  at  one  temperature  can  be 


superimposed  on  that  for  another  by  a  constant  horizontal  displacement 
log  (a^)  along  the  log  (t)  or  log  (  v)  axis,  as  shown  in  Fig.  7. 


In  crystalline  solids  (notably  metals)  the  time-temperature  equivalence 

for  the  rate  of  diffusion,  for  creep,  and  for  other  thermally  activated 

processes,  is  well  known.  The  explanation  is  straight  forward.  Take 

diffusion  as  an  example.  The  frequency  with  which  an  atom  jumps  from  one 

lattice  site  to  the  next  depends  on  the  product  of  two  probabilities:  the 

probability  =  exp  -  (E./RT)  that  the  diffusing  atom  has  enough  energy 

J 

to  clear  the  barrier  (height  E.)  which  separates  the  two  sites;  and  the 

J 

probability  p2  =  exp  -  (E^/RT)  that  the  second  site  is  vacant,  since  if  it 
is  occupied  the  jumping  atom  cannot  enter  it.  The  jump  frequency  is  then: 


v  a  v  exp 
o 


E  E 

-  (J)  exp  -  (j£)  =  ^  exp  -  (gf) 


(7) 


where  Q  =  S  .  +  E  is  the  activation  energy  for  the  process ,  and  v  is 
j  v  o 

the  atomic  vibration  frequency  times  various  entropic  and  geometric  terms 
which  don't  affect  the  arguement.  A  change  of  temperature  from  T  l  to  Tt 
causes  a  change  in  frequency  from  v t  to  v2,  where: 

v  t  Oil 

log  (— )  -  log  aT  =  273-3  -  t7}  (8) 

A  simple  shift  along  the  frequency  or  time  axis  by  log  (a^j  then  brings 
the  response  at  T  k  into  coincidence  with  that  at  T* .  This  particular 
time-temperature  equivalence  is  a  consequence  of  the  operation  of  a  simple 
thermally-activated  process  with  a  single  activation  energy. 

A  spectrum  of  activation  energies  (which  is  inevitable  in  an  amorphous 
system)  does  not  destroy  the  time-temperature  equivalence  though  it  may 
change  its  form.  If  the  width  of  the  spectrum,  aQ,  is  small  (so  that 
AQ  »  RT)  then  it  is  easily  shown  that  the  shift  factor  defined  by  eqn.  (S) 
still  applies,  and  the  relaxation  follows  simple  Arrhenius  kinetics.  The 
3,  y  and  $  relaxations  do,  in  fact,  give  linear  Arrhenius  plots;  those 
for  PMMA  and  PS  are  shown  later.  Sut  the  a  ,  or  glass-rubber  transition 
is  more  complicated.  Experiments  on  liquids  and  on  amorphous  polymers 
(Williams  et  al.,  1955)  are  not  well  described  by  eqn.  (8);  they  suggest. 


instead,  a  shift  factor  (the  "WLF"  shift  factor)  given  by: 


where  C,  and  C,  are  constants  (the  "WLF"  constants,  with  T  as 

g 

reference  temperature).  In  the  limit  T  >>C^  -  T  ,  this  reduces  to  the 

simpler  eqn.  (8)  with  C,  =  Q/ 2.3  RT  .  But  the  values  of  C,  are  such 

g 

that  this  is  never  a  satisfactory  approximation  and  a  new,  broader,  physical 
interpretation  must  be  sought. 


The  quest  had  produced  a  number  of  models  (for  brief  reviews,  see  Ward, 
1971;  Arridge,  1975;  Young,  1981).  From  these  the  following  simplified 

picture  can  be  assembled  and  related  to  the  classical  treatment 
of  thermally  activated  processes  in  crystals  just  described.  Viscous  chain 
motion  is  caused  by  the  stress-biased  diffusion  of  polymer  chains  within  the 
tubes  which  define  their  surroundings ,  as  shown  in  Fig.  8.  This  snake-like 
diffusive  motion,  or  ’’reptation”,  must  involve  the  propagation  of 
compression  or  shear  wave  pulses  or  kinks  along  the  chain:  in  either  case, 
the  projecting  side  groups  of  one  chain  must  move  over  and  past  those  with 
which  they  mesh  in  the  surrounding  tube,  as  shown  in  Fig.  9. 


The  frequency  of  the  unit  step  associated  with  producing  viscous  strain 
depends,  as  before,  on  the  product  of  two  probabilities:  the  probability 
p  j  that  unit  segments  of  the  polymer  chain  have  enough  thermal  energy  to 
jump  over  the  energy  barrier  which  separates  them  from  an  adjacent  position; 
and  the  probability  p2  that  this  second  site  is  "vacant",  i.e.  that  it  has 
sufficient  unoccupied  volume  associated  with  it  to  accommodate  the  jumping 
segment.  The' first  probability,  p  ,  is  given  by  an  Arrhenius  law  like 
eqn.  (7).  The  second  probability,  p  ,  is  more  complicated.  In  crystals, 
unoccupied  or  "free"  volume  is  quantized  as  vacancies  of  fixed  volume.  Free 
volume  exists  in  amorphous  systems  too,  but  it  is  obviously  not  quantized  in 
the  way  it  is  crystals.  When  the  specific  volume  of  an  amorphous  polymer  is 


plotted  against  temperature  (Fig.  10)  there  is  a  change  of  slope  at 


T  . 
g 


The  free  volume  (Cohen  and  Turnbull,  1959;  Haward,  1973)  is  the  difference 


between  the  total  volume,  V,  and  that  occupied  by  the  molecules 


themselves,  V  .  The  occupied  volume  is  that  of  a  dense  (though  disordered) 
packing  of  cylinders,  one  surrounding  each  chain,  with  a  radius  equal  to  the 


Van  der  Waals  radius  plus  that  associated  with  the  local  thermal  vibrations 


of  the  atoms.  What  is  left  over  is  "free"  in  the  sense  that  it 
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redistribute  continuously.  Instead  of  a  thermal -equilibrium  concentration 
of  vacancies,  each  of  fixed  volume,  which  move  around  continuously,  the  free 
volume  exists  (in  thermal  equilibrium,  of  course)  as  a  spectrum  of  void 
volumes  which  continuously  open  up  and  close  again.  The  thermal  expansion 
data  show  that  (unlike  crystals)  the  free  volume  increases  linearly  with 
temperature .  Defining  the  fractional  free  volume  as: 

V  -  V: 

f  m  _ 0 

V  V  (10) 


we  have  that: 


a,  (T  -  T  ) 
f  g 


(11) 


where  f  is  the  fractional  free  volume  at  the  glass  temoerature.  The 
g 

quantity  is  the  free-volume  expansion  coefficient:  it  is  the  difference 

between  the  total  expansion  coefficient  and  that  for  the  occupied  volume. 

The  figure  shows  how  the  value  of  a^.  changes  at  T  ;  above  T  it  is 

large,  below  it  is  much  smaller  (some  authors  take  it  to  be  zero,  but  we 

find  in  Section  3  that  a  small  finite  value  may  be  more  realistic).  The 

probability,  p;f  ,  is  then  the  chance  that,  adjacent  to  a  jumping  segment,  a 

local  free  volume  of  f  or  greater  is  available  (f  is  the  fractional 

c  c 

volume  required  to  accommodate  the  jumping  segment) .  This  probability 
(3ueche,  1953)  is: 


(12) 


Then  the  viscosity  at  temperature  T,  relative  to  that  at  T^  is  given  by 
p  1  (T)  Pj  (f)/pt  (T  )  p2  (fg)>  or: 


n  (T) 

n  (T  ) 

g 


'T 


-  ^r-) }  exp  (3 
g 


(13) 


where  B  =  Af  .  Well  below  T  the  free  volume  is  almost  indeoendent  of 
c  g 

temperature  (Fig.  10);  then  the  first  term  is  dominant  and  we  find 
Arrhenius  behaviour  (though  with  a  spectrum  of  activation  energies).  At  and 
just  above  Tg ,  on  the  other  hand,  f  increases  rapidly  with  temperature; 
and  it  is  probable  that  0,  correspondingly,  decreases;  then  the  second 
term  is  dominant.  experiments  on  viscous  liquids  support  this 
view:  Doolittle  (1951,  1952/  presents  data  wnicn  are  well  described  by 


n  (T  ) 

g 


=  exp  {  3  (—  —  — )  } 


with  the  constant  B  close  to  unity.  Substituting  eqn.  (11)  leads 
immediately  to  the  WLF  equation: 


n  (T) 

lo«v log  run 


-  (B/2.3  f  ) (T  -  T  ) 


f  /a  +  T  -  T 
g  i  g 


with  C  -  B/2.3  f  and  C,  =  f  /a..  One  of  the  many  achievements  of  the 
1  g  g  i 

work  of  Ferry  and  his  co-workers  (Williams  et  al.,  1955;  Ferry,  1961)  is  the 

demonstration  that  C .  and  C,  are  universal  constants  at  and  above  T  . 

1  *  g 

The  same  constants  describe  polymers  also,  implying  that  f  and  a  ,  too, 

g  * 

are  universal  constants.  It  is  helpful  to  note  their  magnitudes:  taking 

B  =  1,  the  fractional  free  volume  at  the  glass  temperature,  f^,  is  found 

to  be  0.025;  and  the  free-volume  expansion  coefficient  above  T  ,  a  -  is 

S  f 

4.8  x  10  /K.  Both  values  are  physically  plausible.  But,  as  before,  the 

simple  rheological  model  of  Fig.  4(a)  leads  to  a  modulus  which  decreases 
too  steeply  with  temperature .  A  much  better  fit  to  the  data  is  given  by  the 
distributions  of  Fig.  4b,  introducing  (as  before)  one  new  parameter:  the 
standard  deviation  <tr  =  ACi/Cj,  and  using  it  in  the  same  way  that 
AQ/Q  was  used  earlier. 


The  Rubbery  Modulus 


Above  T  ,  the  modulus  of  linear,  amorphous  polymers  often  show  a 
plateau  at  around  1  MPa.  This  is  close  to  the  modulus  of  weakly 
cross-linked  rubbers  and  arises,  as  far  as  is  known,  in  a  similar  way.  In 
the  rubbery  state,  the  weak  intermolecular  (Van  der  Waals)  bonds  have 
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largely  melted-  The  long  chain-like  molecules  can  GAB2,  Disc  9 
assume  a  variety  of  configurations  in  response  to  the  thermal  vibration 
which  causes  a  "micro-Brownian"  motion  of  the  units  of  the  chain.  The  most 
probable  configurations  are  those  which  maximise  the  entropy  of  the  system. 
Eut  mere  is  a  constraint;  the  molecular  -cnains  curl 
and  twist  around  each  other;  in  places  they  form  mechanical  entanglements 
which  behave  very  like  chemical  cross-links  (Fig.  11).  When  strained,  the 
chains  tend  to  order  (by  lining  up)  and  the  entropy  decreases;  they  are 
prevented  from  wriggling  back  to  a  disordered  state  by  the  entanglements . 

The  rubbery  modulus  is  related  to  this  change  of  entropy  by  an  equation 
which  is  the  analog,  for  entropy-induced  elasticity,  of  eqn.  (2): 

„  -  d*S  (  e ) 

S  -  T  ~d? - 

where  S  (e)  is  the  entropy  per  unit  volume  as  a  function  of  tensile  strain 


Standard  texts  (Treloar,  1958;  Ward,  1971;  Young,  1981)  summarise  the 
calculation  of  the  rubbery  modulus  from  this  entropy-dominated  model.  For 
small  strains  the  result  is: 


E 


3  g  (T)  RT 
M 


e 


(16) 


where  p  (T)  is  the  density  of  the  rubbery  polymer,-  R  is  the  gas 
constant,  T  the  temperature  and  Ma  is  the  average  molecular  weight 
between  cross-links  or  (here)  entanglement  points. 


The  statistical  mechanics  of  entanglements  in  linear  polymers  is,  at 
present,  beyond  the  scope  of  rubber  theorists;  M  cannot  be  calculated  and 
is  derived  instead  from  data  for  E  by  using  eqn.  (16);  values  are  listed 
in  Table  2-  This  means  that  the  entropy  model  of  rubber  elasticity  has  not, 
strictly  speaking,  been  verified  for  linear  polymers,  although  it  is 
difficult  to  visualise  another  cause  for  the  rubbery  behaviour;  and  the 
values  of  derived  from  the  equation  are  physically  sensible. 

At  first  sight,  eqn.  (16)  suggests  that  the  modulus  in  the  rubbery 

regime  should  increase  with  increasing  temperature,  and  for  many 

cross-linked  rubbers,  it  does.  But  the  density  decreases  with  increasing 

temperature.  If  the  density  at  some  reference  temperature  (say,  T  )  is 

3 

p  (T  )  and  the  volumetric  coefficient  of  thermal  expansion  is  a  ,  then 
3  v 


the  density  at  some  higher  temperature  T  is: 


p  (T) 


p  (T  ) 

_ -8 - 

1  *  a  (T  -  T  ) 
v  g 


(17) 


This  change  in  density  (which  is  included  in  our  calculations)  cancels,  to 

some  extent,  the  explicit  temperature  dependence  of  eqn.  (16). 

Further,  it  is  likely  that  the  weaker  entanglements  unravel  as  the 

temperature  is  increased  so  that  M  ,  too,  increases  with  temperature. 

e 

The  combined  influence  of  this  and  of  thermal  expansion  leads  to  the  plateau 
(usually  with  a  slight  negative  slope)  seen  in  experiments. 


Rubbery  and  Newtonian  Viscous  Flow 

At  high  temperatures  (T  ^  1.2  T  )  the  Van  der  Waals  bonds  melt 
completely  and  even  the  entanglement  points  slip.  This  is  the  regime  in 
which  thermoplastics  are  moulded;  the  polymer  behaves  like  a  viscous  liquid. 
For  a  given  polymer  system  (such  as  PMMA)  the  viscosity,  n  ,  depends  on 
the  molecular  weight,  and  the  temperature  T.  The  time-and- 

temperature  dependent  modulus  (at  constant  stress)  is  simply: 


Two  regimes  of  flow  have  been  identified.  Immediately  after  the  rubber 
plateau,  the  polymer  flows  in  a  way  which  has  the  same  time- temperature 
equivalence  as  that  of  the  visco-elastic  regime  (eqn.  9)  implying  that  it  is 
controlled  by  free  volume.  This  is  generally  called  "rubbery  flow".  Then 
the  viscosity  is  given  by: 


(T  -  T  ) 

n  “  no  SXP  ’  Cj  +  T  -  T 
2  £ 


(19) 


At  higher  temperatures  and  low  shear  rates,  flow  becomes  Netwonian  viscous, 
and  follows  an  Arrhenius  law  with  a  narrow  spectrum  of  activation  energies, 
so  that: 

Q 

n  =  exp  (20)  .  | 
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The  value  of  depends  on  the  molecular  weight,  In ^ ,  in  one  of  two 

ways.  Below  a  critical  molecular  weight,  M  ,  the  viscosity  is 

cr 

proportional  to  molecular  weight  so  that: 


where 


n  is  the  viscosity  of  a  melt  with  molecular  weight  M 
cr  *  cr 


above  M 


the  dependence  changes  such  that: 


"w  3.4 

ncr  (r-) 

cr 


Commercial  polymeric  systems  are  in  this  second  regime;  and  in  the  range  of 
temperatures  covered  by  the  diagrams  shown  in  the  next  section,  rubbery  flow 
dominates . 

The  understanding  of  this  behaviour  is  still  incomplete .  Progress  in 

physical  modelling  has  been  made  by  Ooi  and  Edwards  ( 1978a, b,c,  1979)  who 

extended  work  of  Rouse  (1953)  and  of  De  Gennes  (1971)  to  describe  the 

process  of  polymer-chain  reptat:.on.  The  polymer  chains  behave  like  tangled, 

flexible  strings,  each  enclosed  in  a  tube  made  up  of  the  surrounding  mass  of 

polymer  (Fig.  8) .  Extension  or  shear  of  the  polymer  requires  the  diffusive 

motion  of  the  string  in  its  tube.  In  the  Rouse  (1953)  approach,  the 

3rownian  motion  is  resisted  by  a  Stokes  friction; .equivalently,  compressive 

or  tensile  kinks  are  formed  on  the  chain  and  diffuse  along  its  length 

(Fig.  9):  the  passage  of  one  kink  along  the  entire  chain  displaces  it  by  the 

kink-strength ,  b.  If  the  chains  are  short  and  straight,  as  in  Fig.  9,  the 

time  required  for  a  kink  to  diffuse  along  a  chain  is  simply  proportional  to 

the  molecular  weight.  But  if  the  chains  are  long  and  tangled,  the 

(stress-based)  random  walk,  or  reptation,  of  a  chain  in  its  tube  is  more 

complicated.  Doi  and  Edwards  ( 1978a, b,c,)  and  Doi  (1983)  show  that  the 

— n 

viscosity  should  then  vary  as  in  fair  agreement  with  the  observations 

described  'oy  eqn.  (22). 

Data  for  the  quantities  n  ,  M  ,  M  are  listed  in  Table  2. 

o  w’  cr 


Decomposition 


If  the  polymer  gets  too  hot,  the  thermal  energy  exceeds  the  cohesive 
energy  of  some  part  of  the  molecular  chain  causing  depolymerisation  or 
degradation.  Some  (like  PMMA)  decompose  into  monomer  units;  others  (PE, 
for  instance)  randomly  degrade  into  many  products.  It  is  commercially 
important  that  no  decomposition  takes  place  during  high  temperature 
moulding,  so  a  maximum  safe  working  temperature  is  specified  for  each 
polymer;  typically,  it  is  about  1.5  (Table  2).  We  have  truncated  the 
deformation  diagrams  at  this  temperature. 


3.  DATA  FOR  PMMA  AND  PS 


In  the  last  section,  we  discussed  constitutive  laws  for  each  of  the 

four  main  deformation  regimes.  We  now  fit  data  to  these  laws,  extracting 

values  of  the  parameters  (like  E  ,  a  ,  Q  ,  AQ,  T  ,  etc.)  which  best 

o  m  m  g 

describe  the  data.  The  final  choice  of  parameters  has  already  been 
presented  in  Table  2.  But  the  method  used  to  obtain  them  is  an  important 
part  of  the  process  for  constructing  the  maps  shown  later  •  The  data 
are  drawn  from  many  sources  ,  referenced  in  the  text 

and  on  the  figures. 


Thermal  and  Structural  Data 

Glass  temperatures,  measured  calorimetrically  and  by  dilatometry,  are 
listed  in  Table  3.  We  have  selected  the  values  shown  in  Table  2  as  the  most 
reliable.  The  thermal  expansion  coefficients  (Table  2)  are  from  Yannas  and 
Luise  (1982),  Williams  et  al.  (1965)  and  Van  Xrevelen  (1976).  The  molecular 
weights  given  in  Table  2  are  typical  values  for  commercial  PMMA  and  PS, 
though,  in  practice,  they  vary  widely  depending  on  the  supplier  and  the 
grade . 

The  Glassy  Modulus  and  the  Secondary  Relaxations 

The  modulus  at  absolute  zero,  E  ,  and  its  temperature  dependence, 

o 

a  ,  were  obtained  by  extrapolation  from  above  d  K,  using  the  data  of 
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Sinnott  (1959,  1962).  The  values  obtained  in  this  way  for  PMMA  are  close  to 
the  average  for  all  the  available  data;  and  those  for  PS  are  supported  by 
data  from  Bondi  (1968),  (Table  4). 


The  mean  activation  energy  3  for  each  of  the  secondary  relaxa- 

tn 

tions  i3  found  by  plotting  the  log  of  the  time  t  (or  frequency  v  ) 
for  the  centre  of  the  transition  (or  of  the  damping  peak)  against  1/T.  It 
is  found  (e.g.  McCrum,  Read  and  Williams,  1967)  that  the  data  follow  a 
simple  Arrhenius  law: 


v 


1 

t 


«  exp- 


(23) 


so  that  the  plot  gives  a  stright  line  with  slope  -  0  /2.3  R.  For  PMMA 

(Fig.  12a)  the  3-relaxation  has  a  mean  activation  of  121  kJ  mol”1  ,  in 

agreement  with  the  measurements  of  Iwayanagi  and  Hideshima  (1953  a,b)  and  of 

Sato  et  al.  (1954).  The  activation  energy  for  the  3 -relaxation  is 

13  kJ  mol”1  .  For  PS  (Fig.  12b)  the  8-relaxation  has  a  mean  activation 

energy  of  132  kJ  mol"1  .  Only  two  data  points  (Crissman  et  al.,  1965; 

Sinnott,  1962)  are  available  for  the  5-relaxation  in  PS  so  Q  was 

m 

calculated  from  these  directly. 

When  the  transitions  are  well  separated  (as  they  are  for  PS)  the  drop 
in  modulus  AE  can  be  measured  directly.  When  this  is  not  so  (as  for 
PMMA),  A  E  for  each  transition  is  adjusted  to  give  the  best  fit  for  the 
overall  drop  in  modulus.  Finally  the  spread  of  activation  energies, 

AQ/Om  ,  is  chosen  by  trial  and  error  to  match  the  breadth  of  the 
transition.  The  parameters  which  fit  various  groups  of  data  are  listed  in 
Table  5  and  abstracted  in  Table  2.  The  predictions  of  the  rheological 
model,  using  these  parameters,  are  compared  with  raw  data  in  Figs.  13  (a), 

(b)  and  (c).  In  practice,  the  5 -transition  in  PMMA  and  the  y-transition  in 
PS  have  a  very  small  effect  on  the  moduli  (AE  <0.1  GPa)  and  they  are 
ignored  in  constructing  the  diagrams. 
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The  G las a— Rubber  Transition  and  the  Rubbery  Modulus 

The  glass-rubber  transition,  too,  is  modelled  by  the  rheological 
model  shown  in  Fig.  4(b).  The  modulus  drop  AS 
is  the  difference  between  the  glassy  modulus  (reduced  by 
the  modulus  drops  associated  with  the  secondary  transitions)  and  the  rubbery 
modulus.  The  mean  value  of  Ct  and  its  standard  deviation 

were  chosen  as  described  earlier.  For  PMMA,  a  standard 
deviation  of  0.055  gives  a  good  fit  to  data  (Fig.  14a);  for  PS,  the  value 
0.08  gives  a  good  fit  (Fig.  14b). 

The  rubbery  modulus  is  related  to  the  molecular  weight  between 

entanglements  M^,  by  eqn.  (16),  The  most  recent  measurements  for 

commercial  polymers  are  those  of  Seitz  (1979)  and  it  is  his  data  that  we 

have  used  to  select  Mg  for  PS.  But  for  PMMA  a  range  of  values  for  the 

rubbery  modulus  and  density  have  been  reported.  Results  derived  from  these 

are  listed  in  Table  6;  the  value  we  have  selected  for  M  is  about  half  the 

e 

value  given  by  Seitz.  The  density  p  (T  )  at  the  glass  temperature  was  calculated 
oy  extrapolating  data  for  p  frcm  the  measurement  tenperature  to  T  . 


Figure  15  shows  the  experimental  shift  factors,  log  a^,  for  PMMA  and 

PS.  Above  T  they  are  well  fitted  by  the  WLF  equation  (eqn.  9)  with 

values  of  C,  and  C,  which  are  consistent  with  values  of  f  and  a„ 

S  ^ 

listed  in  Table  7.  The  figures  show  that  at  T  there  is  a  discontinuity 

in  slope  of  the  shift  factor.  The  expansion  data  given  in  Table  7  suggest 

that  is  smaller  below  T  ,  but  that  it  is  not  zero.  We  have  therefore 

^  g 

used  the  WLF  equation  below  T  ,  with  a  new  value  of  C2  ,  calculated  from 
the  data  in  Table  7.  The  final  choice  of  Cj  and  C2  are  listed  in 
Table  2.  The  viscosity  in  the  viscoelastic  transition  has  been  calculated 
as  described  in  Section  2,  using  these  parameters  in  the  WLF  equation. 


Rubber  (Glass)  -  Viscous  Transition  and  Decomposition 


•  .'O 

.*  •  *  o 


Data  for  viscous  flow  are  summarised  in  Table  8.  The  modulus  is 
calculated  by  using  eqn.  (19)  with  the  data  listed  in  Table  2,  including  a 
spread  of  the  constant  C!  Decomposition  data  are  summarised  in  Table  9. 


L.-l 
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4.  DEFORMATION-MECHANISM  DIAGRAMS 

Construction  and  Features  of  the  Diagrams 

The  previous  sections  have  described  the  broad  regimes  of 
deformation  behaviour  for  linear  amorphous  polymers,  between  absolute  zero 
and  the  decomposition  temperature.  Physical  models  exist  for  each 
deformation  mode,  leading  to  constitutive  laws  which  describe  the  time 
and  temperature  dependent  modulus  E;  the  values  of  the  parameters  which 
enter  the  laws  are  known,  at  least  approximately  (Table  2). 

We  now  ask:  over  what  range  of  time  and  temperature  is  a  given 
mechanism  dominant?  And  when  do  the  changes  of  mechanism  appear?  For  each 
mechanism,  the  constitutive  equation  takes  the  form: 

E  *  f  (t,  T,  material  parameters) 

' or  ( 24 ) 

E  «  f  (v,  T,  material  parameters) 

where  t  is  the  loading  time  (or  v  is  the  frequency)  and  T  is  the 
absolute  temperature.  The  dominant  mechanism  for  a  given  T  and  t  is  the 
one  which  leads  to  the  lowest  value  of  E.  Mechanism  changes  take  place 
along  the  lines  obtained  by  equating  pairs  of  the  constitutive  laws.* 


Figures  16  to  19  show  deformation  mechanism  diagrams  for  PMMA  and  PS, 
constructed  from  the  constitutive  laws  of  Section  2,  using  the  parameters  of 
Table  2,  combined  in  the  way  we  have  just  described  .  and  detailed  in  the  Appendix. 
They  shew  the  time-  and  tenperature-depandent  modulus,  E,  for  a  nurfeer  of  different  loading 
conditions,  as  a  function  of  temperature  and  frequency  of  leading. 


Figures  16  and  17  show  the  dynamic  modulus  of  PMMA  over  the  full  ranges 
of  temperature  and  frequency  that  can  be  realised  in  practice.  In  Fig.  16, 

♦The  response  of  the  rheological  model  for  the  relaxations  is  evaluated 
numerically,  summing  the  Gaussian  distribution  of  activation  energies  over 
3i0  on  either  side  of  Q  . 


the  axes  are  E  and  T;  the  contours  are  lines  of  constant  frequency.  The 
four  regions  are  roughly  distinguished  by  the  range  of  modulus  which 
characterises  them,  as  suggested  by  Tobolsky  (1960),  though  strictly  it  is  a 
field  of  values  of  E,  T  and  v  ,  not  merely  of  E,  which  characterises  a 
mechanism.  In  Fig.  17  the  axes  sure  E  and  \>  (plotted  inversely  to  make 
the  diagram  comparable  with  diagrams  which  have  time  as  abscissa:  see 
Gilbert,  1984).  Here,  too,  the  four  mechanism-regimes  are  distinguished  by 
a  range  of  values  of  E.  Figures  18  and  19  show  similar  diagrams  for  PS. 
They  closely  resemble  those  for  PMMA. 


Influence  of  Polymer  Chemistry  and  Molecular  Weight 


It  seems  probable  that  these  diagrams  for  PMMA  and  PS  are  broadly 
typical  of  those  for  linear  amorphous  polymers.  By  normalising  the 


temperature  scale  by  T  ,  the  lowest-order  effect  of  differing  polymer 


g 


chemistry  are  removed. 


At  any  more  detailed  level,  of  course,  there  are  effects  of  polymer 
chemistry  and  molecular  weight.  In  the  glassy  regime  the  secondary 
transitions  are  determined  by  the  nature  of  the  side-groups,  since  these 
influence  the  packing  of  molecules  in  the  amorphous  state  and  so  define  the 
width  of  the  transitions.  (  The  experiments  of  Fujimo  et  al.  (1961),  on 
co-polymers  of  PMA  and  PMMA,  for  example,  show  that  the  larger  the 
sidegroups,  the  more  difficult  it  is  to  pack  the  chains,  and  the  broader  is 
the  glass-rubber  transition.)  The  extent  of  the  rubbery  regime,  too, 
depends  on  the  molecular  weight  of  the  polymer.  Reducing  the  molecular 
weight  shifts  the  contours  in  the  viscous  flow  regimes  to  the  left,  and 
reduces  the  extent  of  the  rubbery  plateau  or  removes  it  altogether. 
McLoughlin  and  Tobolsky  (1952),  for  instance,  find  that  PMMA  with  a 
molecular  weight  of  3600  kg/mol  shows  a  pronounced  rubbery  plateau,  while 
that  with  a  molecular  weight  of  150  kg/mol  shows  none.  Other  studies  of 
the  extent  of  the  rubbery  regime  can  be  understood  in  these  terms 
(Vinogradov  et  al  ,  1971).  And  the  viscosity  of  the  melt  regime,  too, 
depends  on  molecular  weight  (Van  Krevelin,  197«),  though  the  differences 

scale,  broadly,  as  T  . 

3 

We  expect,  therefore,  that  linear  amorphous  will  be  described  by 
diagrams  like  those  of  Figs.  15  to  19,  with  small  differences  caused  by 
molecular  weight  and  architecture. 


5.  SUMMARY  AND  CONCLUSIONS 
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Maps  can  be  constructed  which  summarise  the  time-and-temperature 
dependent  modulus  of  amorphous  polymers,  E  (t,T),  for  a  wide  range  of 
temperatures  and  times,  under  various  loading  conditions.  Several  separate 
mechanisms  are  involved:  bond  stretching,  constrained  molecular  movement  and 
larger  scale  molecular  sliding,  rubbery  behaviour  constrained  by 
entanglements  and  true  viscous  flow.  Each  mechanism  can  be  modelled  (though 
with  differing  levels  of  physical  realism)  to  give  constitutive  equations 
which  describe  how  the  modulus  E  (t,T)  depends  on  temperature,  time  or 
frequency  of  loading ,  and  on  material  parameters  which  characterise  the 
chemistry  and  molecular  architecture  of  the  polymer.  We  have  assembled 
material  parameters  for  PMMA  and  for  PS  and  used  them,  with  the  constitutive 
equations,  to  construct  the  maps  shown  as  Figs.  IS  to  19. 

It  is  probable  that  the  maps  shown  here  are  broadly  typical  of  linear 
amorphous  polymers  (though  maps  for  others  can  readily  be  constructed  using 
the  same  method) .  The  same  approach  can  be  adapted  to  describe  amorphous 
thermosets  and  rubbers ,  and,  with  further  changes ,  to  commercially 
significant  semicrystalline  polymers  such  as  Nylon  and  PE. 

The  maps  shown  here  describe  small-strain  behaviour.  The  next  step,  we 
believe,  is  to  develop  a  parallel  approach  for  large  strain  deformation 
which  is  capable  of  including  cold  drawing,  shear  banding,  twinning  and 
crazing. 
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PMMA 

378a) , 

378b  ^  , 

377C } , 

377d),  388®  \ 

378f } , 

373g  ^ 

PS 

373h) , 

373l) , 

389J )  , 

373g)  , 

a)  Loshaek  (  1955);  b) 
Kovacs  (  1969  );  (d)  Fox 

f)  Van  Krevelen  (1976); 
i)  Plazek  (  1965  );  (j) 


Rogers  and  Mandelkern  ( 1957  );  c)  Wittmann  and 
et  al.  (1958) :  e)  Berry  and  Fox  (1967) ; 

g)  Fox  et  al.  (1976);  h)  Fox  and  Flory  ( 1950) 

Schmeider  and  Wolf  (1953). 


TABLE  4:  YOUNG  ’  S  MODULUS  AND  TEMPERATURE  DEPENDENCE* 


Polymer 

Modulus,  E 

GPa  ° 

a 

m 

PMMA 

8.57a) 

0.28(a 

10.51b) 

7 . 70C  ] 

- 

PS  „ 

6.21c) d) 

0.28d) 

5.80b) 

6.25C) 

Eo  ds  ^oun8  s  Modulus  extrapolated  to  0  K;  am  is  the  normalised 
temperature  dependence  (eqn.  3). 

Sinno^  (1959 )  Extrapoiated  from  4  K  and  eclated  from  G  assuming 


b) 


Yannas  and  Luise 
(d)  are  included 


(1982)  Numerical  average  from  other  sources  (c)  and 
in  the  PS  average. 


c)  Bondi  (1968). 


d)  Sinnott^l962)  Extrapolated  from  d  K  and  calculated  from  G 


assuming 


TABLE  5;  SECONDARY  TRANSITION  DATA 


Polymer 


Transition 


( kJmol" 


‘.Hi 

{ no  a  \ 


AQ  -1 
( kj  mol  ) 


3 

- - - 

iCSjpSp: 

2831} 

r 

I3.0r) 

100°* 

<5 

<20.9s) 

<4.2n^ 

a 

125. 7^J 
136. 3i} 

300° ) 

y 

33.5j) 

132J) 

5 

1 

9.4k) 

36k) 

13  the  mean  activation  energy,  T*  the  temperature 
characterising  the  dartping  peak  at  a  frequency  of  1  !fe ,  i  E  the  modblus 
drop,  associated  with  the  transition  and  AQ  the  spread  in  Q  required 
to  fit  the  experimental  data. 

a)  Deut3ch  et  al  1954 

b)  Iwayanagi  and  Hideshima  1953a,  1953b 

c)  Heijboer  1956 

d)  Sato  et  al  1954 

e)  McCrum  and  Morris  1964 

f)  Sinnott  i960 

g)  Johnson  and  Radon  1972 

h)  Yano  and  Wada  1971 

i) . Connor  1970 

J)  Illers  and  Jenckel  1959 

k)  Calculated  fran  s)  and  t) 

J)  Illers  and  Jenckel  1959 

k)  Calculated  frcm  3)  and  t) 

l)  Heijboer  1965 

m)  Powles  and  Mansfield  1962 

n)  Sinnott  1959 

0)  Illers  and  Jenckel  1953 

p)  Best  fit  to  data  this  study 

q)  Crissman  et  al  1964,  1965 

r)  Schmieder  and  Wolf  1953 

s)  Crissman  et  al  1965 

t)  Sinnott  1962 

u)  Bondi  1966 

w)  Hendus  et  al  1959 

x)  Estimated  value  using  the  method  described  in  the  text 


TABLE  6:  MOLECULAR  ENGANTLEMENT  WEIGHT,  DENSITY  AND  EXPANSION  CQEFr  ICIjiNTS 


Pol ymer 


(kg/ mol) 


9-15b 

10.0 


( 103  kg/m3  ) 


(x-1) 


“v  p(Tg} 

(K-1)  (103  kg/m3  ) 


1 .183“ 
•1.195* 
1.1881 
1.150* 
1 .21  6.h, 


2.7*10- 


6.2*10-'* 
6.6*10-'  ' 


19'1c)  1.0H-1.P55J)  I.SJI  4>5iJ  1 

13.1  1.05fi\  2.8  ;  7.1L) 

1 .057e  . 

1 .052-1 .065* J 

Mg  is  the  molecular  weight  befween  entanglements,  calculated 
from  eqn.  (16),  p  is  the  density  at  the  temperature  listed 


in  the  reference, 


is  the  bulk  exDansion  coefficient  below 


is  that  above  T  ;  p  (T  )  is  the  mean  density 


extrapolated  to  T^. 


a)  Seitz  (1979) 

b)  Masuda  et  al  (1970) 

c)  Onogi  et  al  (1970) 

d)  Fox  et  al  (1958) 

e)  Kolb  and  Izard  (1949)  at  273  K 

f)  Gall  and  McCrum  (1961)  at  298  K 

g)  Rogers  and  Mandelkern  (1957)  at  T_ 

h)  Fujino  et  al  (1962)  at  298  K 

i)  Brady  and  Yeh  (1971) 

j)  Natta  (1955) 

k)  Williams  and  Cleereman  (1952) 

l)  van  Krevelen  (1976) 

m)  Calculated  from  data  presented 
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ia  the  preexponential  viscosity  in  eqn.  (19) 
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TABLE  9:  DEPOLYMERISATION  TEMPERATURES 
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APPENDIX:  COMPUTATION  OF  THE  DEFORMATION  DIAGRAMS 


GABApp,  Disc  3 


The  modulus  E  is  calculated  as  a  function  of  temperature  T  and 
loading  frequency  y  ("dynamic  loading")  as  described  below.  Calculations 
for  modulus  as  a  function  of  temperature  and  time  ("stress  relaxation") 
follow  the  same  method.  The  text  describes  models  for  the  mechanisms 

of  deformation  which  determine  the  modulus,  E.  At  the  highest 

temperatures,  S  is  determined  by  viscous  flow.  The 

modulus  E  _,  is  then  determined  by  the  faster 
vH 

of  the  two  flow  mechanisms  (  eqns  .  19  and  20).  As  the  temperature  is 
decreased,  E  _  increases  until  it  exceeds  the  rubbery  modulus  E_  (eqn. 

VR  ft 

16);  then  this,  not  viscous  flow,  determines  the  modulus:  E  is  set  equal' 
to  the  lesser  of  E  and  E_.  Further  decrease  in  temperature  introduces 

Vft  ft 

the  glass-rubber  transition.  Let  the  increase  in  modulus  associated  with 

this  transition  be  AE  (defined  below) .  Further  drop  in  temperature 

a 

freezes  out  the  3,  y  and  <S  relaxations;  let  the  increase  in  modulus 

caused  by  these  be  E„,  i£  and  a£, .  Then  the  modulus,  at  temperature 

S  Y  o 

T  and  frequency  w  is: 

E  (T,w)  *  Least  of  (E_  and  E  _)  >  aE  +  AEa  +.  &E  +.  &E  (Al) 

ft  Vft  a  °  Y  o 

To  proceed  further,  we  require  expressions  for  the  individual  terms. 
Consider  the  <5  relaxation  as  an  example.  It  is  described  by  the  array  of 
n  Maxwell  elements  shown  in  Fig.  4(b).  For  a  single  element: 

.  a  <j 

€  a  -  +  — 

E  n 

Then  for  the  parallel  array  of  Fig.  4(b),  the  following  constitutive  ec^jaticns  hold. 


Far  constant  strain: 


5  E.  t 

E  ( T ,  t)  a  I  5  E;exp  (-■——) 
i=l 


for  constant  strain  rate: 


E  (T,t)  =  l  -J-(l  - 


5  E  t 

exp  (  -  — - )  ) 

a  n  . 

i 


and  for  dynamic  loading  (so  that  £  »  sin  wt) : 


GLASSY  PLATEAU 


* 


Schematic  of  viscous  diffusive  motion,  or  reptation,  of  a  polymer 
chain  in  the  tube  defined  by  its  immediate  surroundings. 


The  diffusion  of  compressive  kinks  along  a  polymer  chain,  leading 
to  relative  motion  of  chains . 


Heijboer  1965 

Rahaman  and  Scanlan  1981 


~  Schmieder  and  Wolf  1953 
Z 2  Illers  and  Jenckel  1959 
a  Angelo  et  al  1965 
£  Keskkula  and  Turley  1969 
3  Cigna  1970 


Dynamic  Loading 


NORMALISED  TEMPERATURE  T/Tg 


"  '  i  ,  *  3L  t*.  A 

'  * 

t>  Schmieder  and  Wolf  1953 

Z  liters  and  Jenckel  1959 

Z  Angelo  et  ai  1965 

PS 

♦  Bondi  1968 

Dynamic  Loading 

Keskkula  and  Turley  1969 

£  Cigna  1970 

NORMALISED  TEMPERATURE  T7Tg 


Comparison  of  data  from  various  sources  with  the  model: 

(b)  ?S  at  1  Hz  through  T 
(c)  PS  at  1  Hz  at  low  temperature. 


o 

Fujino  et  al  1961 

WLF  Adaption 

V 

WLF  Equation 

V 

PMMA 

TEMPERATURE  /  K. 


Q  Fujita  and  Ninomiya  1957 
71  Grandine  and  Ferry  1953 


WLF  Adaption 


WLF  Equation 

/ 


TEMPERATURE  /  K. 


'he  WLF  shift  factors  calculated  using  the  method  of  the  tex 
:ompared  with  experiments  (a)  for  PMMA,  and  (b)  for  ?S. 


MODULUS  (MPa)  g  g.  MODULUS  (MPa) 


TEMPERATURE  (K) 


NORMALISED  TEMPERATURE  T/Tg 


diagram  for  PMMA  under  dynamic  loading  conditions , 

T/T  as  axes. 

2 


LOG  (FREQUENCY) /Hz 


A  deformation  diagram  for  PMMA  under  dynamic  loading  conditions, 
with  2  and  frequency  v  as  axes.  (The  leg  v  scale  is 
reversed  to  allow  comparison  with  Figs.  20  and  22.] 


g.  IT. 


TEMPERATURE  (  K  ) 


NORMALISED  TEMPERATURE  T/Jq  ■ 


Fig.  18. 


A  deformation  diagram  for  PS  under  dynamic  loading  conditions, 
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A  deformation  diagram  for  PS  under  dynamic  loading  conditions 
with  E  and  frequency  v  as  axes.  (The  log  v  scale  is 
reversed  to  allow  comparison  with  Figs.  21  and  23.) 
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